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The processes of muons (tau) and charged pions pairs production at electron- 

positron annihilation with 0(a) radiative corrections are considered. The calculation 

results are presented assuming the energies of final particles (cms implied) to be not 

' far significantly from threshold production. The invariant mass distributions for the 

\D ' 

C*~) ■ muon (tau) and pion pairs are obtained both for the initial and final state radiation. 

(N , 

*Q \ Some analytical calculations are illustrated numerically. The pions were assumed to 

be point-like objects and scalar QED was applied for calculation. The QED radiative 
corrections related to the final state radiation, additional to the well known Coulomb 

^ ■ factor, are treated near threshold region exactly. 

<D 
43 



I. INTRODUCTION 

The current precision of the evaluation of hadron's contribution to anomalous magnetic 
moment of muon is mainly driven by the systematic error of the cross section of pion pairs 
production at the region where the total cms energy of pair does not exceed threshold value 
significantly. Therefore, the lowest order radiative corrections (RC) and effects due to the 
Coulomb interaction in the final state become essential. 

For the purposes of comparisons with experimental data the cases with hard additional 
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photon must be calculated in frames of PT. It is a weak point of approaches based on 
dimensional regularization methods where the separation of soft and hard photon emission 
can not be arranged. One of the motivations of this paper is to calculate these contributions 
in frames of traditional QED approach with assigning to photon small mass and calculate the 



virtual, soft and hard photon contributions separately. In papers |1| the main characteristics 
of photon emission at annihilation e + e~ to pair of charged particles was investigated. In 
papers published in 1983,1985 the spectra and total cross sections were obtained, but 
the calculation method was too complicated. 

Below using the invariant integration method we repeat in part those calculations and 
obtain the explicit expressions for the spectra distributions on the effective mass of pair 
and the corresponding contributions to the total cross sections due to photon radiation by 
initial or final particles. We do not consider the interference of these amplitudes assuming 
the experimental set-up to be charge-blind. In this case the interference contribution to the 
total cross section is zero. 

In section ITT ] IIHI we consider the final state and initial state radiation of virtual and real 
photons in muon pair production process. In section IIVI similar calculations for the charged 
pion pair production (assuming pion to be point-like object) are done. The results presented 
in sections ITTVlIVI are in agreement with ones obtained in previous papers , but have 

the form more convenient for different applications. Some of them concerning initial state 
radiation are new ones. In section IIVI we also discuss some possibilities of experimental 
separation of contribution of initial and final state radiation. In section the discussion of 
accuracy of results obtained is given. Whenever possible, the analytical results are used as 
a cross-check with ultra relativistic limit. 



II. FINAL STATE RADIATION (FSR) IN MUON PAIR PRODUCTION 

As well as we are interested in muon effective mass spectrum let us put the cross section 
in form: 

spins 

The summed over spin states matrix element squared can be put in form (for notations see 
Fig. EJ): 



e (p_) 




Fig. 1: Final state radiation corrections to e + e" — > process. 
£ |M| 2 = -(47ra) 2 ivT^, S = (p + +p_) 2 , 
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where x± — %kq±, p_ + p + = q = q_ + q + + k, q± = M 2 , p± = m 2 and k 2 = 0. m, M 
- are electron and muon masses correspondingly. Introducing the energy fractions of final 
particles we have: 

2qq± 2qk 

y± — ; v— — , v + u + + v- = 2, 



J dT =J 
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4a p 2 - v 



1 — v V 1 — v ' s 
Due to gauge invariance of tensor T^ u we can write down the following: 



J dTT lllf = \(g lllf -^-\ I dVT?. 



T 

Further simplification follows from gauge invariance of initial leptons tensor U^q^ = 0. 
Simple calculation gives 



A B r, f 

+ + C + (z/+ -> v-) 



'1-V+) 2 '1-1/4 



(3) 
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A = --(3-/3 2 )(l-/3 2 ), C = -2; 
B = -(3- (3 2 ){l + f3 2 ) - 2(3 - (3 2 ) + 2v. 

Integration on the muon energy fraction can be performed using the expressions: 

1 1 



;i-i/+) a ' 



R[y); In tttt' v ™\ v ) 



Distribution on the invariant mass square of muons m 



2 



1 -i?(z/) 



+ 9- 



(4) 



s(l — z/) for the 



case when the energy of hard photon exceeds some value u > y/sA/2, A 1 has a form 
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Fig. 2: Distribution on v for FSR of muons, i.e. the value z = (2a 3 /3s) 1 {da^p SR /dv) (see ©) is 

shown. 
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Contribution to the total cross section can be obtained by performing the integration on 
invariant muon mass. We use the set of integrals: 

/3 2 
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In 



1 + R(u) 
1 - R(u) 
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$(/3) = Li 2 (1 - p) - Li 2 (1 + (3) - Li 2 
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The result is 
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where ctb(s) = 27ra: 2 /3(3 — /3 2 )/(3s) is the cross section of muon pair production in Born 
approximation. In the ultra relativistic limit we have 
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A 4 M 

where = ln(s/M 2 ), £ 2 = tt 2 /6. This result differs from one given in jsj]. The contribution 
of soft real photons emission with energy cu = \Jk 2 + A 2 < y/sA/2, where A - is "photon 
mass", is given by: 



a FSR = a B{S) 
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and performing the standard calculations can be written in form 



a FSR = — a B{s 
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The correction of virtual photon emission include the Dirac and Pauli formfactors of muon. 
It has a form |j| : 

2a l + (3\ M 
a FSR = — o-b(s)[(1 ^— Lp)\n—- 



1 + 



l + (3 2 
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The sum of contributions from virtual and soft real photons reads to be: 
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In ultra relativistic limit we have: 
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v+s 
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The total sum of contributions from virtual, soft and hard real photons does not contain 
photon mass A and the separation parameter A: 
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(17) 



The quantity ApsR. (P) a g rees with the result obtained in |8( and presented in Fig. |3] as a 

function of (3. This correction in ultra relativistic limit tends to the value 3/8. 

47ra; 2 2a 3 a 3 
3s n 8 s 
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e T e — >/j t /x 7| 
FSR 



Cancellation of "large" logarithms = ln(s/M 2 ) is the consequence of Kinoshita-Lee- 
Nauenberg theorem 
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Fig. 3: The dependence of quantity Apjt, (/?) as a function of [3 for FSR of muons. See formula 



(|16l) and it's asymptotic behavior. 
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e + {v+) v + (q+) 

Fig. 4: Initial state radiation corrections to e + e~ — ► process. 
III. INITIAL STATE RADIATION (ISR) IN MUON PAIR PRODUCTION 

Matrix element of the process of muon pair production with hard photon radiated from 
initial state has a form (for notations see Fig. 0J): 



M, 



ISR 



(47TQ! 



,3/2 



s l - i/ 



■v(ph 



p_ - k + m -p + + k + m 
Q e («0 + e ( k ) Q 
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(18) 



— 2fcp_ —2kp + 
with Qr? = u{q-)^ ri v{q + ) is the muon current. Using the gauge condition for muon current 
q v Qn = 0, q = q + + g_ = p + + p_ — k we have: 

E / ^)*^§r^ -«+-*-)= j> - q -f) . (w) 



2tts 



3-/3 2 



7/ 



with notations given above. Using these relations, the calculation of the summed upon spin 

states of matrix element squared is straightforward. Performing the angular integrations 
i 

r i 4 m 2 



dc[ 
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we obtain the distribution on the muons invariant mass (see Fig. |3J): 
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Fig. 5: Distribution of muon pairs as a function on v for ISR. Vertical axis represents the quantity 

y = [(4o 3 /3s)(Z e — l)]^ 1 (da , } SR /di/) (see (|2Tj) h horizontal axis - the fraction of radiated photon 
energy v. 



d °ISR 

dv 



(21) 



3si/(l - z/) 21 

Further integration on the photon energy fraction v can be performed using the set of 
integrals given above and two additional ones: 

/? 2 
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As a result we obtain the cross section due to radiation of hard photon from ISR: 



<TlSR = —<7B(s){l e -l) 
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The contribution to the cross section taking into account the virtual and soft photons to the 
initial state is given by: 



a 



s+v 
ISR 



ix 4 
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Let us note that the spectral distribution on invariant mass of final system have a form 
consistent with renormalization group prescriptions, namely one can recognize the kernel of 
evolution equation contribution (see (J2*T|) . (|2~3Jl ). 
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Now, we can collect all the terms mentioned above and to write out the expression for 
the total cross section due to ISR: 

2a 



a 



s+v+h 



isr ~ T <tb(8)A$ 8 & ((3), 



AC f 



V-W + I3\ 4 2(3 

#3 L/3 ~3 +21n TT^ 



The dependence of this quantity on muon's velocity (3 is shown in Fig. El In ultra relativistic 



3, 

+ 4 h 



(24) 

1 + 6- (25) 




Fig. 6: Distribution of muon pairs on for ISR. See formula (|25|1 for the quantity A^ 5 ^ (/?). 



limit we have: 



s+v+h I 
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which is in agreement with 



DD 



3s -[ 2 Wf 2 lfl 12 /e + 6 ' 24 J ' 
Leading term ~ IJ^ is in agreement with the result 



(26) 



a. 



The total cross section contains the so called double- logarithmical terms (~ Z e Z M ), which 
already contradict the renormalization group predictions (single- logarithmic). 



IV. THE FINAL STATE RADIATION IN PION PAIR PRODUCTION 

It is worth to remind that the total cross section er(e + e - — > 7r + 7r~) with 0(a) corrections 
is required in many subjects of particle physics. Particularly it is required to determine with 
a better accuracy the precision of the evaluation of vacuum polarization effects in photon 
propagator. The other well known application is the calculation of the hadronic contribution 
to the anomalous magnetic moment of muon a^ adr jlfjj ]: 

„ r = i (°Y J ds M£>, R(s) = 4l^lif) (27) 

3 Vtt/ 7 s V ; <r e+e ~^ + ^(s) V ; 

4M| 
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A contribution to this integral coming from high energy region can be calculated within 
QCD framework, while for the low energy range the experimental values R(s) have to be 
taken as an input A numerical evaluation of this integral in relative unities gives the 
value ~ 70 ppm. 

The goal of the new experiment at BNL (E969) is to measure the anomalous magnetic 
moment of muon with the relative accuracy of about ~ 0.25 ppm and to improve the pre- 
vious result |l2j by a factor of two. It follows that the value a l ^ adr should be calculated 
as precisely as possible. In this context the required theoretical precision of the cross sec- 
tions with radiative corrections (RC) as well as the calculation accuracy of the vacuum 
polarization effects should be not worse than ~ 0.2% as it follows from the estimation: 
70 ppmxO.2% ~0.14 ppm. This short observation shows why high precision calculation of 
the hadronic cross sections are extremely important. 



A. Final state radiation 



As well as it was done for the muons, the contributions with one photon radiation in the 
final state can be divided into three separate parts: virtual, soft and hard. The expression 
for the virtual photon emission from final state can be found in |13J and is given by 
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(28) 



Here Lp,\,/3 were defined above, (3 is a pion velocity in cm. frame, 77 (s) = 
(ira 2 p 3 ) j (3s) | F^(s) | 2 is the cross section production of charged pion pair in the Born approx- 
imation, F n (s) - pion strong interaction formfactor. The cross section is due to emission of 
soft photon when its energy does not exceed Ae is given by: 



'FSR 
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Ae < e 
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The sum of the contributions from virtual and soft photons can be presented in convenient 
way as: 

'1 + P 2 
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where 
b(s) 
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Calculations similar to ones given above for muons FSR lead to the pion pair invariant mass 
distribution (m% n = s(l — u), see FigEj): 
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Contribution to the total cross section can be obtained performing the integration on 
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Fig. 7: The pion invariant mass distribution on v for FSR. The vertical and horizontal axes represent 

the value z = [(2a 3 /3s)]~ 1 (da FSR /dv) (see 1(51)1 ) and fraction of photon energy, respectively. 



invariant pion pair mass. In agreement with jj, [lj] the relevant contribution has a form: 
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with $(/?) defined in (jHJ). Now we can write down the complete expression for the total cross 
section: 
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Fig. 8: The dependence of quantity {(3) on (5 for FSR of pions. See formula and it's 

asymptotic behavior. 

with the same expression for F((3) as in muon case (fTTj) . The factor Ajt^ represents 
the correction to the Born cross section caused by final state radiation. In low /3 limit 
Aj.£jj (/?) ~ it 2 /A(3, which is the manifestation of Coulomb interaction of pions. There is 
exactly the same behavior for A^/^ {(3) (see (fbo|) ). It's well known that in limit a/ (3 > 1 
the perturbative analysis is not valid. The relevant modifications of formulae will be given 
in conclusion. 

In ultra relativistic limit we have A^^, (J3 — > 1) = 3/2. One can see again, that all "large" 



logarithms cancelled out in accordance with Kinos 



lita-Lee-Nauenberg theorem. Expression 



for Afcb (P) coincide with one obtained in y, |l5|. It is worth noticing that in papers % 
[hi Q the quantity A^t^ (j3) was presented without separator A between soft and hard 
photons. But for some applications it can be useful to have these two parts separately. 

In order to check experimentally the validity of point-like pions assumption, which is used 
in the paper, it's necessary to separate out the FSR events. Unfortunately we should notice 
that ISR events 10 times more probable then the FSR ones. Nevertheless there are at least 
two ways to select FSR events and to suppress the ISR background. 

Firstly we may consider the region of p-meson peak left slope, i.e. ^/s < 770 MeV. In that 
case the resonance returning mechanism does not take place and the ratio of FSR events 
increases. 

Second way is to throw out the events with pions acollinearity bigger than some predefined 
angle, for instance 0.25 rad [lfil ]. 
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Figure M shows the result of modelling of value uisr+fsr/^isr with application of both 
FSR separation methods described above. The different curves correspond to different energy 
thresholds of emitted photons (u > 10 — 170 MeV). One can see that the energy range from 
720 to 780 MeV is preferable for our purpose - if photon energy exceeds 150 MeV then the 
ratio visr+fsr/visr is about 5, this means that the relative admixture of ISR events is 
about 20 % only. It is worth to notice that the form of spectrum at high photon energies is 




Fig. 9: The ratio of the cross sections with ISR+FSR divided on the cross section with 
ISR as a function of energy in cm. frame. The different smooth curves represent this 
ratio vs threshold photon energy (in MeV) to be detected. 



just the subject of interest. The comparison of the simulated spectrum with experimental 
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one can elucidate the discussed problem. 



B. Initial state radiation 



Let us consider now the initial state radiation (ISR) effects in pion pair production. 
Performing the calculations similar to the case of muon pair production we have: 

a 3 1 + (1 - u) 2 



d °ISR 



-(J, - 1)(P 2 - u)J^^\F n (s(l - u))\\ (35) 



dv 3s (1 — v) 2 v v c " ' V 1 — v 

where q 2 = (q + + q-) 2 = s[l — u). The calculation results are shown in Fig. QUI with = 1. 
Using integrals presented above we can obtain the following expression for the cross section 




0.2 0.4 0.6 0.8 

Fig. 10: The distribution pion pairs as a function on v for ISR. The vertical axis represents the 

quantity y = (a 3 /3s(l e — \))~ l da e Is e R ^ 77 7 /dv (see ([35))). horizontal axis - fraction of radiated 
photon energy. 



with hard photon radiation: 

2a 3 (3 3 
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3s (ie _ 1)iln ^ +21n (_^_ 
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1 l-3[3 2 + A(3 3 

H —n Lp 



(36) 



2(3 2 A/3 3 

where l e = \n(s/m 2 ). Here we had assumed the pions to be point-like, i.e F n = 1. The sum 
of the contributions of virtual and soft photon emission has a form: 
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ISR 



(s) ( (h - 1) In A + -L - 1 + £ 2 

7T 4 



The total cross section accounted for initial state radiation can be presented as: 

2a 3 {3 3 
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1 1 - 3/? 2 + 4/3 3 



+ -/ e -l+6- (39) 



1 + /3 3 2/? 2 4/3 3 
Quantity Af s ^ ((3) as a function of (3 is shown in Fig. [TTJ In ultra relativistic limit in 




Fig. 11: The dependence of the quantity AjJ$ (/3) (see on for ISR. 



point-like approximation for pions we have: 



e^e — >tt 7r 7 
a ISR+FSR 



2a 3 / 1 1 3 1 

/3->i 3s I 2 2 2 6 



(40) 



where = ln(s/M 2 ) 



V. ACCURACY ESTIMATION 



The theoretical uncertainties of the cross sections with 0{a) corrections given above are 
defined by the unaccounted higher order corrections and they are estimated to be at ~ 0.2% 
level. Below the main sources of uncertainties which were omitted in the current formulae 
are listed: 

• Weak interactions not considered here arising from replacement of virtual photon 
Green function by Z-boson one. It results in 



da — > da 



1 + 



Ml) 'Mf 



(41) 



which for ^fs < lOGeV is of order or smaller 0.1% in charge-blind experimental setup, 
when we can omit the 7 — Z interference contribution. 
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Here we systematically omit the terms of order (m/M^) 2 compared to 1 

O gj) < 0.1%. (42) 

The higher orders contributions (not considered here) can be separated by two classes. 
First class, leading by large logarithm l e = ln(s/m 2 ), is connected with ISR: 

da -> da [1 + O {{a/ixf l 2 e , (a/vr) 2 l e )] , (43) 
O ((a/TrR 2 ) ~ 0.2%, O ((a/n) 2 l e ) ~ 0.01%. 



These kind of contributions can be taken into account by structure function approach 
as it was done in 



3- 

• Second class is the higher orders contributions connected with FSR which give 

da l + of^l^j , O^l^j -0.05%. (44) 

In ultra relativistic limit lp —>■ ln(s/M 2 ) they as well can be taken into account by 
structure function method. 

Considering the uncertainty sources mentioned above as independent, we can conclude that 
the total systematic error of the cross sections with 0(a) RC is less than 0.22 %. However 
we remind that taking into account of higher order contributions connected with ISR using 
structure function approach [3] allows one to decrease the total error down to level 0.05 %. 
^fsr an d ^fsR are drawn in figures and one can see that corrections to the Born cross 
sections (2a/7r)A can reach several percents near threshold. 



VI. CONCLUSION 



One of possible applications of formulae given above - to be used for normalization 
purposes at MC simulation. Our results can be used also for improvement of the calculation 
accuracy of vacuum polarization effects in the virtual photon propagator at low energies 
not far significantly from threshold production. This calculation, in one's turn, is required 
to improve the precision of the theoretical prediction for anomalous magnetic moment of 
muon. 
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The expressions for the cross sections of r + r~ and K + K~ production are similarly to 
that for muons and pions. The muon and kaon masses as well as pion form factor should be 
replaced in the above expressions by the tau and kaon ones, respectively. The cross section 
being multiplied by the exact Coulomb factor will interpolate the energy dependence of the 
cross section from the threshold production to the relativistic region. 

We do not consider C-odd interference in real and virtual photons emission - it gives zero 
contribution to the total cross section. As well we do not consider effects of virtual photon 
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3- 



polarization operator insertion, it can be found in literature 

^fsr an d ^isr are drawn in figures and one can see that corrections to Born cross 
sections (2a /n) A can reach several percents near threshold. 

In regions where (3 ~ a formulae must be modified Taking into account, that 
A(*>(/3) ~ vr 2 /4/3, (3 -> 0, we must replace 

l + fA«0<^(l + |(A«0W)-g)), W 

where f(z) = z/(l — e~ z ) is the Sommerfeld-Sakharoy factor, z = (ira/ ' 0). In region there 
,3 C a the formulae must be modified according to 
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